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ABSTRACT 


A  method  is  presented  for  unifying  the  analysis  of  various  wave 
properties  of  a  plasma  in  a  magnetic  field.  An  expression  is  derived 
for  any  microscopic  perturbation  quantity  as  an  integral  of  an  expecta¬ 
tion  value  with  respect  to  the  equilibrium  distribution.  This  yields 
permittivity  and  conductivity  tensors,  and  hence  the  dispersion  relation, 
or  spatial  and  temporal  decay  or  growth  rates,  for  any  specified  velocity 
distribution.  In  the  plane  wave  case,  the  averaging  is  eliminated  and 
the  calculation  significantly  simplified  by  transformation  to  "inverse 
velocity  space,"  so  that  singular  integrals  do  not  appear  and  phenomena 
such  as  Landau  damping  become  evident.  Quasistatic  and  exact  dispersion 
relations  are  derived  for  cold,  Maxwellian,  resonance,  and  drifting 
distributions.  Collisions  are  accounted  for  as  a  viscous  drag  force 
along  the  orbits.  Generalizations  to  other  external  force  fields  are 


discussed. 
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Table  I:  Distribution*  in  Inverse  \  tloeity  Space 
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Figure  1  -  D  taped  helical  trajectorv  of  ‘'avrrage  part  tela 

in  aagnetopl&ana  . . . .  ** 

Figure  2  -  The  oblivion  fa 'tor  due  to  collision* .  22 
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end  spatial  dtoijf  or  growth  rates  are  obtainable  directly.  The  quasi- 
atatic  approximation,  although  often  convenient,  is  not  necessary  to 
the  analysis  and  both  quaslstatlc  and  exact  results  are  presented. 

Singular  integrals,  requiring  careful  specification  of  the  inte¬ 
gration  con t jure  in  the  complex  plane,  are  avoided,  leaving  straight¬ 
forward  quadratures  that  are  either  standard  integrals  or  types  readily 
handled  by  reef* ters.  Landau  and  cyclotron  damping  effects  appear  lu 
a  e tiers I,  enforced  manner,  without  complex  contour  integration  or 
calcelatiee  of  residues  at  ele«4tarlt!es. 

Pinna  ease  or  teener  analysis  i#  performed  la  a  particularly 
einple  mnmrr  If  trenofmcmetien  of  iho  velocity  distribution  to  "inverse 
wM>Hf  epsM/'  «n  •**«*  significant  algobralc  simplification  la 

Its  effect*  ft  drift  or  booms  to  tome  readily  deierrinable. 
fh»  sddoso'i*  !♦  rmilf  gamsfolisod  to  syetoem  etiereelly  forced  other- 

1004*  nop  laeorpmeeie  parametric  effects. 


tie  field. 


2.  SOLUTION  OF  BOLTZMANN  EQUATION 


I  ““ 

The  Boltzmann  equation  governing  the  velocity  distribution  f  of 
a  constituent  of  the  plasma  may  be  expressed  as 


at 


+ 


V  •  Vf  ♦  a  • 


a t 
av 


(2.1) 


Although  it  is  standard6  to  separate  the  collisional  rate  of  change  as 
a  forcing  term,  the  effect  of  collisions  will  here  be  incorporated  in 
the  acceleration,  a  ,  experienced  by  the  plasma  constituent,  as 
detailed  in  Section  4. 

The  nonlinear  equation  will  be  simplified  at  the  outset  by  linear¬ 
izing  about  an  equilibrium  distribution  fQ(v)  .  Thus,  the  distribution 
function  will  be 


f  =  fo  +  fl(£»  *  (2.2) 

where  fQ  conforms  to  the  externally  imposed  force  field  producing  the 
acceleration  ^  in  the  constituent  particles,  and  the  perturbation 
fl  is  associated  with  the  internally  induced  fields  producing  particle 
acceleration  .  With  some  additional  effort,  the  more  general  case 

of  a  spatially  varying  equilibrium  distribution  fQ(r,  v)  could  be 
treated  as  well. 

Linearization  prescribes 

afo 

io  *  %  =  0  (2.3) 

a  v 
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as  a  precondition  for  equilibrium,  and 


+  v  .  rjt  i 


+  -2<> 


(2.4) 


as  the  equation  governing  the  perturbation  of  the  distribution.  How  the 
collision  term  is  to  enter  these  two  equations  is  discussed  in  Section 
4. 

In  these  equations,  a^jr,  ir,  t)  is  a  prescribed  external  excita¬ 
tion  and  the  equilibrium  distribution  fQ(v)  is  presumed  known.  The 
first-order  acceleration  a^r,  v,  t)  is  specified  in  terms  of  the 
r.f.  fields  in  the  plasma,  which  in  turn  depend  on  the  unknown  perturba¬ 
tion  f^  .  A  self-consistent  solution  of  this  equation  is  required  in 
order  to  determine  the  wave  properties  of  the  medium. 

The  solution  to  the  linearized  Boltzmann  equation  is  obtained  by 
integrating  along  the  unperturbed  trajectory  of  a  constituent  particle 
experiencing  the  acceleration  an  .  Provided  that 


dv 


the  equation  states  that 


(2.5) 


(2.6) 
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or  simply  the  total  time  derivative  along  the  trajectory  specified  b; 


a 


(2.5).  The  solution  which  vanishes  at  t  =  -  »  ,  before  the  perturba¬ 

tion  sets  in,  is 

t  if,  /  . 

fl<£,  v,  t)  =  -  J*  -  •  aJr(T),  v(t),  tJ  dT  .  (2.7) 

■oo  Jv  '  ' 

The  integrand  is  evaluated  along  the  orbit  given  by  (2.5).  As  this 
orbit  is  to  pass  through  r  with  velocity  v  at  the  time  t  ,  the 
integral  depends  on  r  and  v  as  well  as  t  .  It  is  implicit  in 
(2.7)  that  the  perturbation  has  a  starting  time.  A  steady-state  analysis 

n 

can  be  reconciled  with  this  by  tacitly  including  a  small  loss  component, 

D 

or  an  adiabatic  switching  factor,  to  guarantee  convergence  of  the 
improper  integral. 


3.  PERTURBATION  QUANTITIES 


The  various  quantities  of  physical  interest  associated  with  the 
plasma  are  expressible  as  expectation  values  of  certain  functions  of 
velocity.  Quantities  that  vanish  at  equilibrium  are  obtainable  by 
averaging  with  just  the  perturbation  f  (r,  v,  t)  ,  typically  yielding 
the  perturbation  quantity  $.(r,  t)  as  the  average  of  a  function  0(v)  : 

X 

o 

$  (r,  t)  =  J*  f  (r,  v,  t)  0 (v)  dv  .  (3.1) 

The  excess  charge  p(r,  t)  is  obtainable  in  this  way  with  0(v)  =  q 

and  the  r.f.  current  density  £(r,  t)  with  0(v)  =  qv  .  More  generally, 
0 (v)  may  represent  a  velocity-dependent  operator  acting  on  the  coordin¬ 
ates. 

In  view  of  (2.7),  any  such  quantity  may  be  calculated  as 

t  dfo 

*1<£.  t)  =  -  J  J  -  •  £1  0 (v)  dT  dv^  .  (3.2) 

-oo  9  V 

This  expression  may  be  simplified  by  invoking  the  divergence  theorem  in 
velocity  space,  along  with  the  condition  that  fQ  -*  0  ,  strongly,  as 

|  v  |  -•  oo  ,  leaving 

*1<£’  t}  =  /V  'oW  fe'[lvFT  '  £i  0(£)J  d£  tiT  •  (3.3) 
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*'o^J  “  •*  *<’»  *»ll.  ik.  wiocl.,  u 

r*«p«ct  to  the  variable  v  •»<») 

U*ln*  angular  bracket.  to  denote  average.  eith  reepect  to  the 
unperturbed  distribution  function,  u 

/  »„<*)  4.  .  .o  ,  «„(«<.»  .  J  »o<,>  ,(,)  d.  ,  <3.4, 

tho  final  roault  la  expreaelble  aa 
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tho  Introduction  of  tho  unperturbed  trajectory  r<v)  ,  *<,)  M  , 
function  of  r  .  v  i 


4.  UNPERTURBED  ORBIT 


Th*  general  result  aay  now  be  specialized  to  the  case  of  a  constant, 
uniform  magnetic  field  as  the  external  force  field.  This  introduces 
s  preferred  direction  in  space,  that  of  ^  ,  to  be  specified  by  the 
unit  vector  b  .  It  is  then  convenient  to  define  a  trio  of  matrices 
associated  with  b  ,  as  follows. 

j|.bb,  -  b  b  ,  X  =  b  x  .  (4.1) 

**•**•»  .£  unlt  »*trixj  ||  and  J_  are  seen  to  be  projection 

operators  and  X  performs  the  cross  product  operation.  A  mutually 

orthogonal,  ideapotent  set  of  matrices  is  formed  by  R  ,  L  ,  || 

where 


R  -  jtl-iX)  ,  L  =  i  (1+  iX)  .  (4.2) 

These  represent  right-  and  left-handed  circular  polarization  operators, 
respectively.  The  spectral  expansion  of  the  operator  X  shows  that 
any  matrix  function  of  X  reduces  to 

f(X)  .  f(i)  R  +  f(-i)  l  +  f (0)  ||  .  (4.3) 

The  utility  of  these  definitions  arises  from  the  fact  that  in  the 
external  magnetic  field  ,  the  constituent  particle  acceleration 


produced  is 


=  (q/m)  v  x  £0  =  -  X  v  ,  (4.4) 

t 

where  coc  =  qBg/m  is  the  signed  cyclotron  frequency.  With 

=  dv/dt  ,  this  equation  describes  the  precession  of  the  velocity 
vector  about  the  magnetic  field. 

One  consequence  of  the  form  of  the  acceleration  imposed  by  the 
external  magnetic  field  is  that  the  condition  (2.3)  on  the  unperturbed 
velocity  distribution  becomes 


!fo 

o 

/O 

<• 

=  -9. 

.  X  •  V  =  D  •  V  X  - 

dv 

Sv 

dcp 

which  precludes  any  velocity-space  azimuthal  (cp )  dependence  of  the 
equilibrium  distribution. 

Although  the  introduction  of  the  orbit  into  (3.5)  is  now  straight¬ 
forward  for  the  collisionless  case,  it  is  desirable  at  this  point  to 
incorporate  the  effects  of  collisions,  in  some  manner.  This  would  at 
least  serve  to  resolve  ambiguities  associated  with  singularities  appear¬ 
ing  in  the  absence  of  collisions.  Some  of  these  arise  from  the  fact 
that  the  matrix  X  is  singular.  To  maintain  tractability,  collisions 
are  here  to  be  included  in  the  simplest  fashion,  in  terms  of  an  equiva¬ 
lent  "collision  frequency"  or  inverse  relaxation  time,  v  .  This  is 

commonly  introduced  In  any  of  various  convenient  or  reasonable  approxi- 
9-11 

mations  to  the  collisional  term  in  (2,1).  It  should  be  noted  that 
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•  simple  relaxation  tara  In  tba  Boltzmann  aquation  is  Inadequate  when 

the  quaaistatic  approximation  is  to  be  used  ab  initio,  due  to  Its 

o 

failure  to  conserve  particles  locally. 

The  artifice  to  be  employe''  here  to  represent  coilislonal  effects 
is  to  ascribe  them  to  morely  a  modification  of  the  unperturbed  orbit. 

The  typical  particle  is  considered  to  be  subject  to  a  viscous  drag 
force,  in  addition  to  the  magnetic  one,  to  the  extent  that  its  velocity 
differs  from  the  mean  flow  velocity  at  equilibrium.  Besides  the  resul¬ 
tant  tractabillty,  this  method  has  the  virtues  of  yielding  results 
consistent  with  the  limiting  case  of  cold  plasma  hydrodynamic  theory, 
as  well  as  consistency  between  quasistatic  and  exact  theory.  Although 
this  approach  neglects  diffusion  in  velocity  space,  the  proposed  change 
in  the  unperturbed  orbit  appears  also  as  an  essential  modification  of 
the  collisionless  case  when  more  careful  account  is  taken  of  collisions 
via  a  Fokker-Planck  model. 

Accordingly,  the  acceleration  in  (4.4),  which  prescribes  the  orbit, 
is  modified  to 


—0 


=  -  O)  X 


0 1  -  v  -  Zo>  • 


(4.6) 


where  v  is  the  effective  collision  frequency,  assumed  constant,  and 
v^  =  (v)  is  the  mean,  d.c.  drift  velocity  of  the  equilibrium  distribu¬ 
tion.  The  acceleration  is  still  linear  in  the  velocity  in  this  model 
and  the  orbit  is  readily  expressed,  from  (2.5),  by 

dr/dT  s  ^(t)  »  dv/dT  =  -(cc^X  +  v)v  +  vVq  ,  (4.7) 
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•luce 


•ubjett  to  tbs  eo*  mot*  £Cf)  *  £  .  £<T>  *  £  •*  T  *  1  * 

V  ii  <o  be  tims-isd#peadeat  Md  Is  necessarily  aligned  »lth  the 

4> 

upttllc  field,  ths  orbit  equation  actually  simplifies  to 


£  <l-.v  “<u,oi*  v)<~“  v  • 


(4.«) 


in  tirai  of  tbs  peculiar,  or  randoc,  velocity  v  -  ^  and  the  nonaingular 

matrix  a)  X  +  v  *  Y  . 
o-  **■» 

The  solution  for  the  trajectory  is 


vCO-Vq  «  e  ~fT't)  (v  -  Vq)  5 


r  +  v^r-t)  +  [1  -  «'~(T_t)]  Y_1(v  -  v^  . 


(4. ft) 


(4.10) 


Equation  (4.9)  expresses  the  damped  precession  of  the  peculiar  velocity 
vector  about  b  ,  until  it  attain#  v  -  at  time  t  ;  eq.  (4.10) 
describes  the  constricted  helical  path  taken  by  the  particle  which 
ax  *ives  at  r  with  velocity  v  at  time  t  ,  as  indicated  in  Fig.  1. 

By  use  of  (4.3),  these  expressions  may  be  rendered  more  explicit, 
for  reference  purposes.  Let  T  «=  t  +  u  ,  £(T)  *  £  +  £  s  £^r)  “  £  * 


£  =  £ Q  +  U(u)  •  (£  -  vQ)  , 


(4.11) 


s  «  vQu  +  S(u)  .  (v  -  v0)  ( 


(4.12) 
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ORBIT 


Figure  1  -  Damped  helical  trajectory  of  "average  particle”  in  magnelopt**** 


rifc)  *  •  -  *  •  (* 


*  **  1  -  *!•*»  •)«  ,  {) 


M.SJI 


and 


-ru%  -i 


»(U)  .  (1  -  .**")  Y 


I  -  • 


■(V*1q  )U 


I  ♦ 


1  -  • 


-(V-Ito  ) u 


1  -  • 


-vu 


v  ♦  lw 


v  -  k±> 


.  (4.14) 


Note  that,  since  vn  ■  I  v 

~0  1  -0  ’ 


XDio  *  «<v>io  (4.1 

tor  any  Mtrlx  fraction.  Al.o,  the  re.ult  s<u)olu  .  -s(-u)  will 
useful  later. 

It  may  be  noted  that,  by  virtue  of  the  time-invariance  of  the 
system,  the  particle  dynamics  are  entirely  expressible  in  terms  of  the 

elapsed  time  u  =  T  -  t  .  Finally,  for  use  in  (3.5),  note  that 

av/&v(T)  =  Ju  . 

Substitution  of  these  results  in  (3.5)  makes  it  applicable  to  the 
magnetized  plasma: 
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•»«<  •»  *S  t 


;4  [• 


*i'tt***  H 


•  * 


^  W  l*l 


**  *  *m  *  *#*  b»*m**»  #|  »  m«i  *?*,** 

»  «»  at 

IMflnllHI  I*  M  tl«lt)  WrffMli*  |l«*»  ito  ft  **«*+<*•* 

HIM!  «f  III*  fMtlllf  f«  |(t)  #p  «•  «<***fl- 

<«**M  «t  <M*f»  «  M4  MH  •  •!  lb*  MfHIIlN  |lM«,  #%)«*  *M*f- 

I*  *.(r«  *,  i)  I*  lb*  itHtN  f«fi«  11*14#.  H*  latur 

^1  4*  *» 

•rt  tbu  ttwiwifwiit  n#U*  *1(11*4  i*  Ik*  #)**•*,  m  that 


*lf£*  *.  tj  •  (*/■)  (tjlr.  t>  .  *#  ijCr,  t>)  .  <4. IT) 

Th*  aourct*  of  *j  ,  »w  ik*  cUrp  ud  currant  danattlas,  obtaln- 
abla  In  turn  froa  (4.16). 
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5.  MACROSCOPIC  OBSERVABLES 


fh*  b»sic  physical  quantities  of  interest  are  the  excess  charge 
an4  curreot  densities  in  the  plasma.  These  nay  be  calculated  by  taking 
0(t)  •  q  and  ^v  ,  respectively,  in  (4.16): 


p(r,  t)  -  oQq  J 


J(r,  t)  *  nQq  J 


r°  <1 .  n 

loo  L 

C4-E 


(to  X+V)u 

•  1' 


(03  X+V)U 

e  °~  •  SjCrf^,  w, 


(5.1) 


(5.2) 


It  may  be  noted  at  this  point  that  the  general  result  of  eq.  (4.16) 
reduces  to  the  following  expressions  in  the  cases  of  harmonic  time 

r  i03t  i 

variation  and  of  plane  waves.  With  s^Or,  v,  t)  =  Re  [s^Oriir)®  J  or 
a  (r  v,  t)  «  Re  [a  (vie1*"*'**-*  ]  ,  and  corresponding  expressions  for 

the  macroscopic  quantities  of  interest, 


(5.3) 


(5.4) 


Wave  properties  of  the  medium  are  most  readily  obtained  from 
quantities  derived  from  the  charge  and  current,  rather  than  from  the 
sources  themselves.  Dispersion  relations  are  expressible  in  terms  of 
susceptibilities  and  spatial  and  temporal  decay  rates  in  terms  of  the 
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power  balance.  To  obtain  this,  Maxwell's  equations,  or  their  equivalent, 
must  be  combined  with  the  above  expressions. 

The  medium  may  be  considered  as  either  dielectric  or  conducting, 
as  convenient.  The  dielectric  susceptibility  K(k,o>)  is  so  defined 
for  plane  waves  that 

P  =  “i€(&  ’  £<£»“*  *  £1  (5.5) 

and  the  normalized  conductivity  C(k,co)  is  defined  by 


£  **  “leoa)  ££»«>  *  Ei 


(5.6) 


The  dispersion  relation  for  wave  propagation  is  conveniently  expressed 
in  terms  of  these  tensors.  Under  the  quasistatic  approximation^,  it 
suffices  to  set 


k 


2 


k  .  K(k,o>)  •  k 


(5.7) 


since  is  considered  to  be  the  gradient  of  a  plane  wave  potential. 

The  exact  plane  wave  dispersion  relation,  however,  is  obtained  by  com¬ 
bining  Maxwell's  equations  into  a  wave  equation,  with  the  current  as 
source .  This  yields 


eig 


k  k 


-  1 


(5.8) 


where  the  notation  eig  A  »  x,  means  that  \  is  an  eigenvalue  of  A 
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The  equivalent  permittivity  and  conductivity  tensors  are  obtainable 
from  (5.4)  by  translating  (4,17)  into  plane  wave  notation.  The  magnetic 
field  may  be  consistently  dropped  under  the  quasistatic  assumption, 
leaving  a^  =  (q/m)E,  ,  but  the  exact  acceleration  is 


"  k  k 

a  (v)  =  (q/m)  (1  -  -  •  v)  _I  +  -  _v  •  E  .  (5.9) 

~~L  ""  L  oo - OO  ""  J  ~ 


2  2 

Substituting  in  (5.4)  and  writing  nQq  /meQ  as  oop  ,  the  square  of 


the  plasma  frequency,  leads  to 


k  •  K(k,oo)  =  ici)  J' 


:4  ■  [■ 


(oo  X+v+ioo)u  -ik. 
o— 

e  e 


~1\ 

y,du  (5. 


10) 


and 


oo 

C(k,oo).Ei  =  i  -EJ* 


•[< 

(JO  -oo  *- 


(O)  X+v+io>)u  -ik 
o— 

e  e 


: »  f  *  *  ^  \ 

ui  —  .  w)_i  +  —  WJ  •  ^ndu 


(5.11) 


Taking  the  divergences  in  velocity  as  indicated  simplifies  these  to 


K(k,oo)  =  oo  2  J  j*(u).e 


(oo  x+v+ico)u  -ik.s 
°~  <e  ~  du  (5.12) 


and 


C(k,<J0)  =  i 


(00  X+v+iOO)u  -ik 
o» 

e  e 


k-s  f  k  k  “|A 

I  (1  -  —  •  w) I  +  -  w  fldu 

L  a  '“w  «)  j/ 


(5.13) 
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In  bounded  systems,  the  quantities  of  interest  are  often  decay 
rates,  spatial  or  temporal,  particularly  when  these  can  be  negative, 
indicating  instability  or  amplifying  capabilities.  From  the  power 
balance  in  terms  of  power  and  energy  densities, 


7 


-J.E 


f 


(5.14) 


there  is  obtained,  by  appropriate  integration  in  space  and  time,  the 
decay  rates  of  energy  in  a  cavity  or  power  in  a  waveguide,  as  follows. 


20L 


W  T 

av 


u  i* 


dV  dt 


(5.15) 


2as  =  pV  I J*  i'l  **  dt 
av  T  A 


(5.16) 


In  terms  of  the  formalism  developed  here,  these  decay  rates  may  be 
calculated  by  taking  0(v)  to  be  qv<E  ,  integrated  in  space  and  time 
The  result  of  using  this  velocity-dependent  operator  in  (4.19)  is 


2a 


o/  r  (<«x+v)u  “V 

=  no  q  *  Le  • 


\du  ,  (5.17) 


where  the  correlation  tensor  is 


R(jB,w,u)  =  if  ^  ^  *i(r+£,  w,  t+u)  £(r,t)  dV  dt  (5.18) 


for  a  cavity,  or 


»(.,  i.  u)  .  jrif  J  ;  ijCr*.  «•  «*»)  £(r,t)  1*  «  (5.U) 


av  t  A 


for  a  waveguide.  Other  aacroacopic  obaervablea,  auch  aa  frequency  ahlfta, 
are  similarly  obtainable  by  appropriate  choice  of  the  aicroacoplc  opera¬ 
tor  0(v) 
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«.  INVKm  VELOCITY  SPACE 


A»  I.  evident  fro.  <5.12),  the  basic  quantity  to  be  averaged  over 
velocity  apace  has  the  fora 


-iki 

e 


-ik-vu  1A  • (v 
•  °  e  “  - 


(6.1) 


.here 


■  ^(u)  •  -k'8(u)  (6.2) 

IKtrine  generally,  therefore,  the  distribution  function  in  "inverse 
velocity  space"  by 


r<A)  .  .  <..3) 

This  is  Just  the  velocity-space,  three -diaens tonal  Fourier  transfora  of 
the  noraallaed  unperturbed  dlatrlbution  function  f0(v)/n0  .  The 
equivalent  peraittivity  tensor  in  (6.12)  is  then  expressible  directly 
in  teras  of  r^A)  .evaluated  **  (*.*).  The  averaging  of  aore  coapll- 

cated  functions  of  velocity  eill  often  be  advantageously  expressible  in 
tvnte  of  f(£)  as  *sll,  as  in 

(*•  *  -■)  *  "l  !  ( vve^**) 


.i!r. 


(6.4) 
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In  particular,  the  various  velocity  moment,  are  easily  obtained  by 
evaluating  the  derivatives  ol  Fft)  at  the  origin  A  =  0  .  the 

Laplaclan  ol  F<A>  at  the  origin  Is  -<»2> 

The  algebraic  form  ol  the  distribution  lunctlon  Is  generally 
simpler  In  Inverse  velocity  space,  because  ol  the  normalisation 
r<0)  =  1  and  the  replacement  ol  convolution,  by  product,,  m  parUcu. 
lar,  a  drilling  distribution  appear,  simply  ..  the  stationary  one 
multiplied  by  the  exponential  lactor  a1*  Mo  :  Table  I  give,  the 
lorn  Ol  various  distribution.  In  both  velocity  spaces.  that  the 

temperature  T  is  taken  in  energy  units. 

In  terms  of  the  distribution  function  in  inverse  velocity  space, 
the  equivalent  permittivity  tensor  becomes 


v  .  2  X+v+ i  (cu-k •  v  )  lu  -iAev 

K(k,a>)  =  &  *  J  £|(u)  e  —  ~Q  J  e  ^  Zc 


F(A)du  ,  (6.5) 


A(u)  _  -k'S(u)  and  “-k.^  iB  the  Doppler-shifted  frequency. 

The  equivalent  conductivity  tensor  becomes,  after  some  manipulation 
of  (5,13), 
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•  m 


TABU  I 


Distributions  in  Inverse  Velocity  Space 


1 .  General : 


2.  Cold: 


f(v)  =  nQ  8  0v) 
3.  Drifting: 


'<i>  -  Vi  -  V 

4,  Maxwellian: 

,  _  mv2 

f(v)  =  n0(m/2jtT)3^2  e  2T 

5.  Resonance : 

n  v 

f  “  ~ 2  ~~2  2?2 
It  (V  +  V^ 


f(v)  =  f(  V  ) 


f(v)  =  f (v  ) 

p 

8.  Isotropic  monoenergetic: 


f  (v)  = 


nos<W  •  V 


4  Jt  v 


9.  Transverse  monoenergetic: 


x  nn8(v  "  vi>&<v-) 

f  (v)  ^  ..  ft,  .  p _ 1  z 

2  it  v 


F(A)  =  J  f(v)e^'~  f (v)d3v 


F(A)  =  1 


FQ\)  =  e  ^  F  (A) 


F(A  )  =  e 


F(A)  =  e 


F(A  )  = 


1  1a2 

7  in  A 


Bin  A  v 


F<4)  =  <J0ft  »  )) 


■W  vi 


F(A )  = 


F(A)  = 


10.  Beam  in  plasma: 


f<*>  ■  Vp<~)  +  nb£b(^o) 


11.  General  averages 


nn"nh 

F<A)  =-VF  p<A>  +  JTe  W 

0  y  0 


(tty)  -  !*&  tCv)  d3v/nQ  =  J  F<A)  *<A>  d\ 
where  *  (A  )  =  (2Jt)'3  J  t  (v)  e"1^- ~  d3v  ^  *  <v)  -  J  t  ^.>  . 
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is  the  contribution  of  the  r.f.  electric  field,  and 


» 


£b  = 


CD 

_P 

2 

CD 


-f  i*  X 


i(^-k-Vo)u 


-^•lo 


-00 


4(- 


-CD  XU 
C— 


SA 


CDgXuX 


X  A  *e  c'”  j  du  (6.8) 


is  due  to  the  r.f.  magnetic  field. 

If  cavity  or  waveguide  modes  are  expanded  in  plane  waves,  similar 
reductions  of  eq.  (5.17)  may  be  achieved.  There  remains  to  substitute 
any  appropriate  distribution  function  F(A)  into  these  expressions  to 
yield  the  dispersion  relations  or  decay  rates  by  straightforward  quadra¬ 
ture.  Singular  integrals  do  not  appear  in  this  formulation. 

More  generally,  the  perturbation  of  the  distribution  function  in 
inverse  velocity  space,  F^^)  ,  is  obtainable  by  setting 

no0(v)  =  e1^  in  (5.4).  Thus,  the  perturbed  distribution  is 

i(CDt-kT) 

F(r,A,t)  =  F(A )  +  F1(A)e  "  ,  (6.9) 


with 


Fi(A) 


i(CDU-^C'£  +  A  ^v)  V 
(w)e  \ 


du 


(6.10) 


The  perturbations  in  charge,  current,  temperature,  etc.  are  obtainable 
by  evaluating  this,  and  its  derivatives,  at  A  *  0  . 


7,  QUASISTATIC  DISPERSION  -  B  =  0 

Under  the  quasistatic  approximation,  the  significant  quantity  is 
the  equivalent  permittivity  tensor,  given  by 


K(k,03) 


0  to  X+v+i  (<u-k  •  v  )  ]u  -iA  •  vn 

J  S(u)  e  e  ~  ~°  F(A)  du 

-00  ^ 


(7.1) 


The  dispersion  relation  is  then  k2  =  k<K(k,co)<k 

The  case  of  BQ  =  0  is  recovered  in  the  limit  <o  =  o  ,  whereupon 

c 

S(u)  =  [(1  -  e  VU)/vJl  ,  or  S(u)  =  u£  in  the  collisionless  case,  and 
the  permittivity  reduces  to  a  scalar:  K(k,oo)  =X(k,o>)  j 


X(k,o)) 


1 


-  e 


-Vu 


V 


[v+i(o>-k*v  )]u  -iA  *v_ 

1^(1  (J 

6  ®  F(A)  du 


(7.2) 


where  A  =  (k/v)(e  -  1)  .  The  dispersion  relation  is  simply 

XOj,^)  =  1  .  Note  that  v^  a  -i  &F(0)/dA 

a)  Cold  plasma:  F(A )  =  1  . 


X(k,<o) 


2  »0  i  _  e-Xu  (V+iu»u  #  • 

%  J  - —  •  = 


u>(<n-iv) 


(7.3) 


A  small  loss  component  has  been  invoked  to  make  the  integral  convergent. 
The  dispei ‘ion  relation  here  represents  merely  s  damped  oscillation  at 


-ilA2 

2  a 


ihe  plasma  frequency. 

b)  Maxwellian  distribution:  F<£)  -  e 
nere,  A2  -  (k/v)2(e"Vu  -  l?2  or  A2  •  k2u2  for  v  .  0 


The 


chilis ionless  c  se  is  the  most  tractable,  reducing  to 


-  >4  - 


X  (k,co) 


_  1  k_T  2 

...2  pO  iaou  2  m  U 
a)  f  ti  e  e  du 


2 

ffl  b  o  „  2 

4-  J  26  e2y6  e~e  d 6 

k  T  -oo 


(S?Ti 

P 

2 

k  T 


dS 

dy 


(7.4) 


where  y  .  (to/kJ^T)1/2  and  E(y)  =  J°  e"5^2^  d5  .  The  integral 

“00 

is  the  error  function,  with  complex  argument.  Landau  damping,12'13  is 
implicit  in  this  result.  It  appears  again  more  explicitly  in  the 
following  case. 

x  -v-,  |A| 

c)  Resonance  distribution:  F(A)  -  e  ,“"1 

Here,  |^|  =  <k/v)<e  ^  "  1)  or  |a|  =  -ku  for  v  =  0  .  The  collision¬ 

less  case  exhibits  damping: 


2  viku 

X^k.co)  «  (Op  J  ue  e  du  =  - £ - -  .  (7.5) 

-*  (<w  -  iv^r 

c  1 1 1  ■  ion  lea  ■  Landau  damping  at  tba  rate  v^k  la  evident  from  the 

dispersion  relation.  The  physical  interpretation  of  this  case  is 

obscure,  however,  as  this  distribution  has  no  well-defined  not**nt* 

lA-v^ 

d)  Beams:  f<A  >  -  e  ~  rQ(A )  . 

k  Staple,  distributed  beam  results  la  only  a  Doppler  shift 

X(h,*>  -Vi,*  *  1'V  * 


.15  . 


as  is 


evident  from  (7.2).  For  .  be-  la  •  .t.tlonery  Pl-~.  h°"'v'r' 


2  ...2  .  m2 


n  _  i  +  and  o>p  in  (7.2)  is 

The0  X0  s  Zb(nb/D0)  and’  in  (7,2>' 


.2)  is  to  be  replaced  by  =  “p  +  ’ 


ii.  --b 

°o  rV  =  “p  FP  +  “b  ®  ~~  F*  * 


(7.7) 


equal,  opposed,  interpenetrating  beams,  -  ®  6,1(1 


F(j\)  =  2u>2  cos  ^*vb  Fb(^)  . 


(7.8) 


• n 


8.  QUASISTATIC  DISPERSION  -  B  y*  0 

In  a  magnetic  field,  £(u)  is  given  by  (4.14).  The  product 
S(u)e^U  appearing  in  (7.1)  reduces  to  -S(-u)  ,  leaving 


0  i  (in-k»v,.)u  ik*S(u)  •  v 

K(k,co)  =  Y_1  J*  (e^U  -  1)  e  ~  e  F(-k*S(u))du  (8.1) 


where  Y  =  v  +  co  X  .  In  the  collisionless  limit,  -S(-u)  reduces  to 
—  o«- 


sin  (0  u  /  1  -  cos  cu  u 

-s(-u)  =  - —  i  +  - — 

V  “c  . 


X.uJ 


(8.2) 


and  the  dispersion  relation,  which  involves 


sin  a)  u 

-k»S(-u)*k  =  -  k  +  u  k„  , 


(8.3) 


reduces  to 


i  -  i  «n 

W  “c  - 


(8.4) 


where 


i  “-e 

o  <A> 

g^  .  J*  sin  0  e  C  F(A)  d 0 


(8.5) 


27- 


and 


0)  „ 

i  —  Q 

U)_ 


g„  =  J°  6  e  C  F(A)  d9  , 

-oo 


with  A  =  -(k/o^)  »  $sin  &1.  -  (1  -  cos©)  X  +  #) 

a)  Cold  plasma:  F(A)  =  1 
From  (8.1), 


K(k,o)>  =  a)2  Y_1  f°  (elU  -  l) 


1o)U  . 

e  du 


0) 


-E  (O)  -  iY)"1 


O)2 

_£ 


0) 


O) 


R  L 


o)-iv+o)  (jD-iv-o)  o)-iv 
•  c  c 


and  the  dispersion  relation  is 


O)  2  (O)  -  iv) 

r  2  i 

P 

.2 

+  kN 

p 

a)[(a)-iV  )2  -  a)  2] 
l_  C  J 

a)(u)-iV)_ 

0  . 


b)  Maxwellian  distribution:  F(A )  =  e 


’*  T  —  A2 
2  m 


In  the  tractable  collisionless  case,  with  9  =  u>  u  , 

C 


R  k 

A2  =  —  •  8(6)  •  S' (6)  •  —  a  2(k  2/o>  2)(1  -  cos0)  +  (k«i2/a)  2)  02 

0)  ~  ~  0)_  1  C  "  C 

c  c 


Let 


X  = 
l 


k  2T 
1 

2  *  1 
a) 
c 


^2l 


X*li  "  n  JT  X  =  X  +  Xfl 
a)czm  " 
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(8.6) 


(8.7) 


(8.8) 


.  (8.9) 


(8.10) 


Then 


F(A) 
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X  cos  B 
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(8.11) 


so  that 

1  y  a2 

o  X  cose  2 

gi  =  e  /  8lne  e  e  e  d6  (8.12) 

"00 

and 

-x  o  X  cose  ‘  2 

=  e  J  e  e  e  A  e  de  ,  (8.13) 

-oo 

where  0  a  o)/^  .  By  Inspection  of  the  Integrands,  it  nay  be  antici¬ 
pated  that  the  parallel  component  of  k  will  introduce  collisionless 

13  7  <- 

damping.  ’  The  special  case  of  perpendicular  propagation  is  free  of 

this  and  is  of  particular  interest.  The  dispersion  relation  is  then 


1  S 


-x  o  me  x  cos e 

~2  e  |  sine  e  e  1  de 


Integrating  by  parts  and  using  the  Fourier  series 
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00 
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yields 


1  k2T 
1  +  2~  =  e 
“p " 


-X 


I0(x)  +  2ft2  Y 


oo 

e  \(X) 


£i  0Z-nS 


(8.16) 


This  dispersion  relation  for  perpendicular  propagation  displays  cyclotron 

harmonics.  It  is  easily  shown  to  agree  with,  but  converge  ‘aster  than, 

the  version  quoted  by  Stix.14  it  is  equivalent  to  that  given  by 
15 


-v. 


.1*1 


Bernstein. 

c)  Resonance  distribution:  F(A)  a  e  '1,~l 

The  collisionless  case  of  perpendicular  propagation  is  readily  handled. 
From  (8.9)  with  k||  =  0  , 


F(A  )  =  e"2M  (sin  0/2 J 


(8.17) 


where  m  =  v^/o^  .  The  dispersion  relation  is 
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which  reduces  to 
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<-1>ni2n(2M) 
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p  V 

L 

“00 

<P+n)2  -  1 
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00 


<-l)“+1<!otirfl  I2n+1<2p) 

..  .1  2  . 

(P  +  n  +  |)  "  1 


.  (8.18) 


This  also  shows  cyclotron  harmonics  but  not,  despite  appearances, 
singular  behavior  when  0  =  n  +  1/2 

d)  Effects  of  drift  and  collisions:  F(/0  =  e~  F(A)  . 

—  0 

A  drifting  distribution  introduces  no  more  than  a  Doppler  shift  into 
(8.1),  as  is  physically  evident  from  a  transformation  to  a  moving 
reference  frame.  If,  however,  there  is  drift  of  one  plasma  component 
relative  to  another,  as  exemplified  by  (7.7),  the  mean  drift  velocity 
dm®rs  from  that  of  either  component.  There  are  then  Introduced  into 
the  Integrand  in  (8.1)  factors  such  as  exp  ik.S(u)*(vrt  -  v.)  besides 
the  Doppler  shift.  A  brief  discussion  of  their  effect  is  in<  order. 

By  (4.15),  the  drift  factors  are  of  the  form 


~k'(vn  "  vh)  vul 

exp  ik-S^).^  -  =  exp  i  - - —  (1  -  e  )  (8.19) 

\)  J 


In  the  collisionless  case,  this  is  exp  ik»(v^  -  v.  )u  ,  which  merely 
reassigns  to  each  drifting  component  its  proper  Doppler  shift.  In  the 
presence  of  collisions,  the  effect  of  this  factor  is  clarified  by  the 
interpretation  of  the  original  integral  (4.16)  as  a  superposition  of  the 
perturbations  of  the  past  (u  <  o)  ,  as  propagated  to  the  present 
(u  ~  0)  .  The  factor  (8.19)  is  an  oscillatory  function  ei^<u>  ,  of 

instantaneous  frequency  cty/du  =  k. (^  -  e_VU  ;  see  Fig.  2.  In  the 
sufficiently  distant  past,  this  is  of  so  rapid  variation  as  to  erase  all 
memory  of  earlier  perturbations,  as  confirmed  by  the  well-known  Riemann- 
Lebesgue  theorem.  This  oblivion-producing  aspect  of  collisions  is  seen 
to  be  incorporated  mathematically  in  the  model  of  collisions  adopted  herein. 

-  '31  - 


m«  exact  dlaperalon  relation,  obtained  frou  the  full  aet  of 
Mexvell'a  equetlone,  ia  given  in  tenaa  of  the  noraeliaed  conductivity 
teneor  C(k,o>)  by  (5.S).  In  thia  connection,  it  ia  uaeful  to  note  the 
following  proportion  of  that  equation,  Firet,  the  electric  field 
ia  to  be  an  eigenvector  correepondlng  to  the  preacilbed  eigenvalue. 
Second, 


eig  (eg  ♦  |)U  •  O  oil  *  ♦  P 


(0.1) 


la  an  identity.  Third,  if  the  conductivity  teneor  abould  take  the  for* 


C  •  o  _I  ♦  P  kk/x 


(t.ti 


the  diape ralon  relation  (S.S)  would  becone 


*“{[(?)  ‘ p]  ?  '  °)  '  (?)  ' 


(S.3) 


or  by  (0.1), 


°  *  (?)  *  [p  ‘  (?)  ]  •“  ? 


(0.4) 
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2 

But  «ig  kk/k  ■  1  for  longitudinal  aodea,  or  taro  for  tranaverae  wxlaa. 
The  diaparalon  relation  then  dcconpoaea  Into 


1  »  ft  ♦  li 


for  longitudinal  aodea,  and 


(8.S) 


1  *  °‘(»)  <»•< 

for  tranaverae  aodea.  The  Matrix  C(k,m)  la  given  generally  by 

(6.6-6) .  If  £  reduce#  to  o  ,  a  actlar,  <6.6.6)  atlll  apply,  but 
•1th  0*g 

In  the  abaence  of  the  external  Magnetic  field,  o>  •  0  and 
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which  i*  of  tb*  fora  (9.2).  Th*  diapersion  relation*  are  therefore 
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for  longitudinal  *ivm,  and 


1  - 


U)  . 

‘  -i  £ 

0)  • 


[i  +  (v^/co)  e] 


de  +  — 


for  transverse  waves.  Explicitly,  these  are,  respectively, 


and 


which  exhibit  Landau  damping. 

*A\£o 

d)  Cold  beam:  F(A_)  =  e 

Here  HE  =  i  v  F  and  both  the  r.f.  electric  and  magnetic 
’  5A_  ~0 

contribute  to  C  ,  yielding  in  the  collisionless  case, 


(9.27) 


(9.28) 


(9.29) 


(9.30) 


fields 


(9.31) 
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10.  EXACT  DISPERSION  -  B  ^  0 


In  the  presence  of  a  magnetic  field,  the  exact  expression  for  the 
normalized  conductivity  tensor  entering  the  dispersion  relation  is,  as 
previously  given,  C(k,cD)  =  C  (k,CD)  +  C.  (k,CD)  where 


O) 


C  (k,co)  =  i 


ro  -VL'Zo 

J  e 


CD 


('♦gi) 


[Y  +  KcD-k-v^u 


du 


(10.1) 


with  Y  =  cu  X  +  v  ,  and 
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CD2  ' 


■“oiu  jf  , 

'  •  5A  XAe 


CD  Xu\  i(CD-k*V„)U 

o-l  v - -0 

/e 


du 


(10.2) 


In  these, 


A(u)  =  -k-S(u)  =  k-Y"1^  ~U  -  lj  .  (10.3) 


The  dispersion  relation  Is 


with  the  electric  field  as  eigenvector. 


(10.4) 
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a)  Cold  plasma:  F(^)  =  1  . 


£<kfa>) 


2 

0)  0  (^  X+v+ia>)u 

i  -E  J  .  O'  du 
0)  “00 


iu>  x) 

O'' 


(10.5) 


The  electric  field  vector  must  satisfy 


c2(o)-iv)  k  k*E  -fa)2  a)  +  (c2k2-  032)(ci)-iv)*l  E  =  c2o>  k-E  iX-k  -  (c2k2-o>2)co  iX. 

(10.6) 

Waves  of  various  linear  and  circular  polarizations  may  be  extracted  from 
this  equation. 

-  -  -  A2 

b)  Maxwellian  distribution:  F(\)  =  e  2  m 
Since  SF/fc^  =  -(T/m)AF  and  X'  =  -X  ,  the  cross  product  in  '^(k,^) 
vanishes,  leaving 
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In  the  cases  of  these  last  two  distributions,  only  special  cases  of 
propagation  and  polarization  parallel  and  perpendicular  to  the  magnetic 
field  can  be  considered  tractable.  It  is  useful  to  note  the  explicit 
expression 


but  the  complexity  of  the  quadratures  precludes  further  development  of 
general  dispersion  relations  here.  Of  considerable  interest,  however, 
are  the  spatial  and  temporal  decay  rates  of  waveguide  and  cavity  modes 
and  the  possibilities  of  amplification  and  oscillation  when  these  can 
become  negative . 
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1  .  A 

R(s,  w,  u)  =  - — -  J  j  ^(r  +  £,  w,  t  +  u)  E(£,  t)  dVdt  ,  (11.2) 

av  TV 

the  integrations  being  over  the  cavity  volume  V  and  a  time  interval 
T  .  In  this  equation,  the  acceleration  a^(r  ,  v  ,  t)  is  the  Lorentz 
one,  given  in  (4,17).  In  (11.1),  s  and  w  are 


-Yu 

w  s  vQ  +  e  — 


,  ,  /i  -Yu,  -1  .  , 

-  Vq)  ,  s  =  vqu  +  (1  -  e  ~  )  Y  •  (v  -  vQ) 


(11.3) 

where  Y  =  o>  X  +  v 
—  o- 

To  evaluate  the  decay  rate  of  a  cavity  mode,  sufficient  accuracy 
is  obtainable  by  using  the  unperturbed  field  pattern  in  these  expressions. 
The  correlation  tensor  may  be  evaluated  for  any  desired  cavity  mode  and 
the  substitution  for  w  and  s  in  (11.3)  then  permits  the  integration 
of  (11.1)  explicitly. 

When  the  cavity  mode  is  readily  decomposed  into  plane  waves,  the 
calculation  of  the  decay  rates  is  particularly  simple,  as  then  the 
distributions  in  inverse  velocity  space  are  directly  utilizable.  The 
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calculations  will  here  be  illustrated  for  the  simple  case  of  the  "cavity" 
consisting  of  infinite  space,  wherein  a  suitable  mode  is  the  plane  wave 
£  =  Re  E^e  .  The  acceleration  is  then 

-2i(£'  Z>  t)  =  Re  a  ["(i  _  r  .  V)  !  +  r  v")  .  e 

1  mL  m  -J  -i  ai.«> 

and  the  correlation  tensor  is  readily  found  to  be 


. 


a  r  k  k  "1 

R(s,  w,  u)  =  ~  Re  I  (1  -  -  .  w)  £  +  -  w  I 

o  L  ^  ~  ~  aj  ““ J 


~i~i 


5r*i 


i(am-k.s) 


(11.5) 


Hence,  by  comparison  with  (5.11), 


2a  =  Re 


-iOD 


E  •  C(k,cn) -E  ” 

“Vi  -V  •*» 

£i*  •  £i 


lm 


£1  •  a>C(k,a)).Ei 

Si*  •  Ii 


(11.6) 


The  results  previously  obtained  for  C(k,a>)  may  therefore  be  used 
directly  to  calculate  decay  rates  in  the  plane  wave  case. 

In  the  absence  of  the  magnetic  field,  the  decay  rate  for  a  cold 
plasma  is,  from  (9.9) 


2a  =  Im 


<n-iv 


o?  v 

P 

2  2 

a>  +  v 

P 


(11.7) 


and  is,  of  course,  due  to  collisions. 
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For  a  Maxwellian  distribution,  the  decay  rate  obtained  from  (9.14) 
is,  for  longitudinal  waves, 


a)2  o 


2a  =  Re  -H  J  (i  -  r2e2) 


202)  e16  e 
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-  5  7  0 

2  2  “ 

70)  cose  e  de 


<4/*fa  .-Y<2- > 


0)  V  2 


which  is  the  collisionless  Landau  damping  decrement,  with  y2 
For  transverse  waves,  similarly, 


(11.8) 

Tk2/ mo)2  . 
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(11.9) 


For  the  collisionless  case  of  a  resonance  distribution,  the  corres¬ 
ponding  decay  rates  are,  from  (9.29,30), 


2a  =  Im 


(05  -  iv^) 


2  2 

2a>  a)  v.k 

P  1 


.  2  2  2  2 
(0)  +  v1  k  ) 


(11.10) 


for  longitudinal  waves,  and 


2a  =  Im  - H_ 


a)  -  iv^k 


mp  v 

cd2  +  Vj2  k2 


(11.11) 


for  transverse  waves. 
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In  the  presence  of  the  magnetic  field,  the  decay  rate  for  a  cold 
plasma  is,  from  (10.5), 

2  **  -1  A 

2a  =  Im  03  e  •  (03  -  iv  -  103  X)  .  e 
p  ■  o~  ~ 

e*  *  R  '  e  e*  *  L  *  e  e*  .  J[  .  e 

(03+03  )2  +  V2  (03-03  )  2  +  V2  032  +  V2 


where  ^e  =£i/)®1|  •  Note  that  iX  ,  R  ,  ,  Jj  are  all  hermitiap 

matrices,  whereas  only  the  antihermitian  part  of  C(k,03)  contributes 
to  the  decay  rata. 

For  a  Maxwellian  distribution,  from  (10.7), 


2a  =  Re  03 


*  x  nO  T 

£  -  J  <£-5*4.3 


(03  X+V+i03)U 
<3- 


1  T  2 

2  m  A 


(11.13) 


This  reduces,  for  example,  in  the  collisionless  case  of  a  parallel- 
propagating  transverse  wave,  to 


(11.14) 
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The  right-  and  left-handed  circularly  polarized  waves  are  seen  to  suffer 
different  Lendau  damping  decrements. 

The  resonance  distribution  case  analogous  to  this  last  one  yields 
similarly 


„  0  (oj  X  +  iau  +  v_k)u 

2  **  p0  e~  1 


2a  =  Re  uj  '  •  J  e 


du  «  o 
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(co-a>  )  +  v.  k 

C  L 


(11.15) 


The  effects  of  a  beam  component  of  a  plasma  in  a  magnetic  field 
are  of  special  interest,  because  of  the  possibility  of  negative  decay 
rates.  An  indication  of  how  this  can  come  about  is  provided  by  the 
following  considerations.  The  inverse  velocity  space  distribution 

Q-Zo 

function  for  a  warm  beam  has  the  form  e  FQ(A)  *  wIlere  Zq  is  the 

beam  velocity  and  F^(A)  gives  its  stationary  distribution.  In  the 

expression  for  C^k.oi)  ,  eq.  (10.1),  F  +  (&F/dA)A  is  hence  replaced  by 
1A  •  Vf.  "r 

e  ~  [fq  +  (S F ^A  M  +  iv^  A  FqJ  .  The  exponential  factor  cancels 
in  the  integrand,  leaving  a  Doppler  shift.  The  last  term,  arising 
from  the  drift,  can  provide  an  imaginary  component  for  Cq (k,<h)  ,  whose 

sign  depends  on  that  of  the  drift  velocity  v^  with  respect  to  the 
direction  of  propagation.  For  example,  the  coniribui.ion  to  C  (k,co)  of 
this  last  term  in  the  case  of  a  beam  wich  a  Maxwellian  distribution 
drifting  along  the  external  magnetic  field  is,  in  the  absence  of  collisions, 


du  , 


(11.16) 
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Q.yc+i(iB-k.v0;ju 


so  that  for  a  parallel -propagating  longitudinal  wave,  for  which  A 
reduces  to  -ku  ,  the  contribution  to  the  decay  rate  is,  from  (11.6), 


2a  _  = 
eO 
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1  V-li-lr/ 

2  (T/m; k2 


(11.17) 


This  can  be  negative,  either  if  vQk  <  0  ,  i.e.  for  upstream  propagation, 

or  if  k.vQ  >  cd  ,  i.e.  if  the  Doppler  shifted  frequency  is  negative. 

This  negative  Landau  damping  decrement  would  have  to  overcome  the  normal 

Landau  damping  provided  by  the  other  terms  in  the  expression  in  order  to 
leave  net  growth. 

The  Fourier  analysis  into  plane  waves  is,  of  course,  not  necessary 
for  thev calculation  of  the  decay  rate  for  a  cavity.  The  correlation 
tensor  R(s,w,u)  is  a  property  of  the  cavity  field  pattern  which  may  be 
calculated  separately  before  it  is  introduced  into  eq.  (11. l). 
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12.  SPATIAL  DECAY  RATES 

The  rate  of  decay  of  power  flow  along  a  waveguide  rr.*y  be  calculated 
in  a  manner  completely  analogous  to  that  of  temporal  decay  rates.  The 
formula  is  the  same, 


2a 


‘  £(£»-w*  u) 


but  the  correlation  tensor  is  now 


(1?.  ,1) 


1 

R(£,  w,  u)  =  - — -  P  J  £-0*8,  w,  t+u)  £(r,t)  dA  dt  , 

av  T  A 

(12.2) 

where  the  integrations  are  over  the  waveguide  cross  section  A  and  a 
time  interval  T  .  In  (12.1),  £  and  w  are  functions  if  v  and  u  , 

as  in  (11.3).  The  correlation  tensor  is  a  property  of  the  waveguide  field 
pattern  and  may  be  calculated  independently,  before  the  substitutions 
for  s  and  w  make  it  a  function  of  £  and  u  . 

Again,  the  decomposition  of  the  field  pattern  into  plane  waves, 
when  convenient,  expedites  the  calculations,  and  again  the  simple 
illustration  of  an  infinite-space  "waveguide"  is  instructive.  For  a 
plane  wave,  the  energy  density  and  power  flow  are  related  by  the  group 
velocity  of  the  wave.  Hence,  the  spatial  decay  rates  for  infinite  plane 
waves  are  obtainable  from  the  temporal  decay  rates  derived  above  by 

2  i 

dividing  in  each  case  by  the  undamped  wave  group  velocity  c  k/o>  . 

For  a  rectangular  waveguide  supporting  the  TE  mode 
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E(r, t)  =  Re  e  sin(k»r)ei propagating  along  £  ,  polarized 

along  e  ,  with  k»p  =  0  and  |kj  =  ir/a  ,  the  correlation  tensor  is 
found  to  be 


R(s,  w,  u)  =  (q/mc2€gP)  Re  Qg  cos  k«s  +  i  h  sin  k«s~]  e  , 

(12.3) 

where 


g  =  e(a>-j3.w)  +  £  e  •  w  ,  h  =  e  k.w  -  k  e.w  .  (12.4) 


Hence,  thi  spatial  decay  rate  is  obtainable  as 


2a  =  (a>^c2|3)  Re  e  .  J°^I  +  v  •  e^+^)U  .  [j  cos (k-s )+ihsin(k 


(12.5) 


—Yu  -1  -Yu 

with  S*/5v  =  e  —  and  3s/dv  s  _Y  (1  -  e  —  )  .  If  the  trigonometric 

terms  were  decomposed  into  exponentials,  the  previously  derived  plane 
wave  results,  in  terms  of  the  distribution  in  inverse  velocity  space, 
could  be  used  directly,  by  appropriate  composition. 
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13.  OTHER  EXTERNAL  FORCES 


The  results  obtained  apply  to  a  plasma  subjected  to  external  forces 
limited  to  a  constant  magnetic  field,  together  with  an  equivalent  fric¬ 
tional  drag  force  to  represent  collisions  phenomenologically.  The 
approach  can  be  readily  generalized  to  allow  more  complicated  force 
fields.  These  might  include  gravity,  d.c.  electric  fields,  time  varying 
pump  fields,  and  incident  waves.  The  latter  two  would  result  in  para¬ 
metric  effects  and  various  wave-wave  interactions. 

To  introduce  any  external  force  field  requires  only  the  solution 
of  the  dynamics  problem  giving  the  unperturbed  orbit.  For  example, 
there  may  be  added  to  the  magnetizing  field  and  viscous  force  a  constant 
electric  bias  field,  or  gravity.  Then 

i 

*o  =  +  +  (q/m)  lo  +  v^o  (13-1) 

and  the  orbit  is  given  by 

v(t)  =  e"~U  ■  v  -  y"1  (e_~U  -  1)  [(q/m)E0  +  V  v J  (13.2) 


and 


_r(T)  =  r  +  Y_1  (1  -  e"~U)  •  -  Y-2(l  -  e"~U)  -  Y_1u  (q/m)K0+Vvo  , 


(13.3) 


where  u  =  r  -  t  .  These  are  to  be  substituted  for  w  ar’d  r  ♦  s  , 


respectively,  in  the  various  equations. 


For  a  time  varying  applied  electric  field  EQ(t),  the  orbit  is 


^(T)  =  Zq  +  e~^U  •  <i-v0)  +  (q/m)  JU  e“-^  E  (t  -  §  )  d?  .  (13.4) 


KJJ  l 

In  particular,  for  a  harmonic  pump  field  E_(t)  =  Ee  0  .  this  reduces 


-Yu  if  icunu  “Yu"l  iOD.t 

I(T)  =  Z0  +  e~  *  (Z~Z0)  +  <2  +  %>*  L®  “  e~  J  (q/®)Ee  .  (13.5) 


1(a>ot“iSn’r) 

If  the  pump  is  an  incident  wave  E(r,t)  =  EQe  f  the 

orbit  is  given  by  the  solution  to 


dr  dr  q  dr  /k^ 

.  v  ~  0 


dT  m 


_+Y  E0+-X  ~XE 


dT  \  a) 


l[V  -  io’-£<T0 


+  v^o  • 


(13.6) 


with  \^(t)  =  dr/dT  and  subject  to  _t(t)  =  r  and  v(x)  =  v  at  T  =  t 
Since  the  entire  approach  is  perturbational,  it  would  not  be  inconsistent 
to  use  orbits  obtained  from  this  equation  by  successive  approximations  or 
linearization. 
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14.  CONCLUSIONS 


A  formalism  has  been  presented  which  unifies  the  analysis  of  the 
various  wave  propertiea  of  a  pi....  subjected  to  some  external  force 
field.,  particularly  magnetization.  The  aim  has  been  to  derive  relations 
of  a  general  nature,  into  which  there  may  be  Introduced  the  appropriate 
descriptions  of  both  the  unperturbed  orbit  in  the  applied  field  and  the 
qulllbrium  velocity  distribution  of  the  plasma  constituents.  Straight¬ 
forward  quadratures  then  yield  the  wave  dispersion  relations,  the  permit¬ 
tivity  or  conductivity  tensors,  the  charge  or  current  distribution.,  and 
absorption  or  growth  rates,  both  temporal  and  spatial,  ..11  ..  the 
perturbed  velocity  distribution. 

In  summary,  the  Boltzmann  equation  has  been  linearized  by  separating 
the  acceleration  into  externally  applied  and  internally  induced  components 
with  collisions  considered  phenomenologically  as  a  viscous  retarding 
force.  The  solution  to  the  linearized  equation  was  taken  beyond  merely 
that  for  the  perturbing  velocity  distribution  to  obtain  directly  an 
expression  for  the  perturbation  of  any  macroscopic  quantity  which  can 
be  calculated  as  an  expectation  value  with  respect  to  the  unperturbed 
velocity  distribution. 

The  effects  of  collisions  have  been  accounted  for  in  a  convenient, 
yet  not  unrealistic,  manner.  With  just  an  effective  relaxation  rate  V 

“  Par“eter.  th'f  »»  considered  as  simply  damping  the  otherwise  helical 
unperturbed  orbits  of  the  constituent  particles.  Collisions  appeared  as 
“  *d<utl°“*1.  effectively  external,  viscous  drag  force  tending  to  relax 
the  velocities  toward  the  mean  flow  velocity.  This  model  i,  mathematically 
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tractable  and  avoids  tbe  paradoxes  associated  with  Indiscriminate  conver¬ 
sions  of  real  frequencies  to  the  complex  form  m  -  lv  .  It  achieves 
consistency  >otween  quasistatic  and  exact  results  and  reduces  properly 
to  the  cold  plasma  fluid  model.  In  the  formulation  Interpreted  as  super¬ 
positions  of  perturbations  along  the  particle  trajectories,  the  collision 
model  introduces  factors  that  tend  to  destroy  the  system's  memory  of 
perturbations  suffered  in  the  distant  past,  as  measured  by  v  ^  , 

There  is  thus  provided  a  simple  mathematical  model  of  physical  collislonal 
effects . 

The  key  quantities  obtainable  in  this  formalism  as  expectation 
values  are  the  permittivity  and  conductivity  tensors,  from  which  the 
dispersion  relation  can  be  extracted  by  combination  with  either  the 
Poisson  equation  or  Maxwell's  equations.  Spatial  and  temporal  decay  or 
growth  rates  were  expressed  in  terms  of  similar  correlation  tensors 

f.  « 

associated  with  the  waveguide  or  cavity  field  pattern. 

In  the  case  of  plane  waves,  it  was  found  that  most  quantities  of 
interest  are  expressible  through  the  expectation  value  of  an  exponential 
function  of  velocity.  This  was  the  basis  for  a  significant  simplifica¬ 
tion  introduced  by  expressing  the  equilibrium  velocity  distribution  in 
inverse  velocity  space;  i.e.  as  a  Fourier  transform  in  velocity.  A 
variety  of  expectation  values  are  then  obtainable  by  simple  differentia¬ 
tion  or  convolutions.  In  addition,  singular  complex  integrations  are 
thereby  avoided  and  phenomena  such  as  Landau  and  cyclotron  damping 
appear  naturally. 

Explicit  results  have  been  presented  under  the  quasistatic 
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approximation  and  without  it,  in  the  absence  and  in  the  presence  of  the 
applied  magnetic  field.  Dispersion  relations  have  been  derived  for 
these  cases  for  cold  plasmas,  for  Maxwellian  plasmas,  for  resonance 
distributions,  and  for  beams  or  drifting  plasmas,  in  some  cases  with 
collisions  included  explicitly. 

Temporal  and  spatial  decay  rates,  particularly  Landau  and  cyclotron 
damping,  have  been  calculated  for  plane  waves  by  simple  integration. 

The  quadratures  necessary  in  less  tractable  cases  have  been  indicated. 
Finally,  the  requirements  for  generalizing  the  theory  to  include  para¬ 
metric  interactions  have  been  presented. 
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